Surprising relations between parametric level correlations and fidelity decay 
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Relations between fidelity decay, cross form-factor (i.e., parametric level correlations), and level 
velocity correlations are found both by deriving a Ward identity in a two-matrix model, and by 
comparing exact results, using supersymmetry techniques, in the framework of random matrix 
theory. A power law decay near Heisenberg time, as a function of the relevant parameter, is shown 
to be at the root of revivals recently discovered for fidelity decay. For cross form-factors the revivals 
are illustrated by a numerical study of a multiply kicked Ising spin chain. 
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Fidelity decay presently attracts considerable atten- 
tion [l|]. It measures the change of quantum dynamics 
of a state under a modification of the Hamiltonian. In 
quantum information, fidelity measures the deviation be- 
tween a mathematical algorithm and its physical imple- 
mentation. From a different point of view, important 
insight into the properties of the underlying systems is 
provided by the studies of correlations between spectra 
of random and/or chaotic Hamiltonians which differ by a 
parameter-dependent perturbation 0]. Since statistical 
properties of fidelity decay in random/chaotic systems 
involve both spectra and eigenfunctions of the original 
and perturbed Hamiltonians, existence of any connec- 
tions between fidelity and purely spectral correlations is 
not a priori obvious. 

Random Matrix Theory (RMT) has been successful in 
describing quantum many-body systems and as model for 
the spectral properties of single particle systems whose 
classical analogue is chaotic |4[. Within RMT fidelity 
was analyzed in linear response approximation [a and 
both fidelity [H @, 0] and parametric correlations Q were 
calculated exactly using the supersymmetry method. An 
unexpected fidelity revival at Heisenberg time was en- 
countered Q within RMT and confirmed in a dynamical 
coupled spin chain model 

Earlier, differential relations between parametric spec- 
tral correlations and parametric density correlations were 
established 1(J HI • By relating the latter to the fidelity 
amplitude via Fourier transform, we show in this letter 
that the existence of these relations opens a crucial in- 
sight into the properties of fidelity decay. By analyzing 
the characteristic features of the parametric correlations 
in the time domain, the cross-form factor, we discover 
a new, simple interpretation of the previously puzzling 
phenomenon of revival [3| . These relations follow directly 
from the basic definitions and symmetries of the underly- 
ing matrix models, being essentially Ward identities. We 
show that they are valid under very general assumptions. 
No explicit (e.g., supersymmetric) calculation is required, 
however they rely on the universality of the parametric 



spectral correlations at the scale of mean level spacing. 
We thus explain the origin of various relations connecting 
spectral and wave-function correlations, and establish a 
unified framework for their analysis and generalizations. 
A relation between fidelity decay and level velocity corre- 
lation function is given. The latter is important from the 
experimental point of view, being used for independent 
access to system parameters. We confirm the general re- 
sults comparing fidelity decay and cross-form factors in 
RMT. We illustrate our analytical results with a numer- 
ical study of a multiply kicked Ising spin chain. 

We consider Hamiltonians modeled by N x N matrices 



H (X) = H±XV/2 



(1) 



where H and V are independently drawn from ensembles 
of the same symmetry. In particular, V is drawn from the 
GOE, the GUE or the GSE ensembles of RMT, labeled 
(3 = 1, 2, 4. The ensemble average over both is indicated 
by angular brackets. It is convenient to fix the variances 
as (HijHki) = D^ 1 (VijVki) where D is the mean level 
spacing of if ± (0) in the energy region of interest. In the 
RMT case, D = ir 2 /N in the center of the spectrum. 
The mean level spacing is then A-independent up to cor- 
rections of order 1/iV. By construction, i? ± (A) is in the 
same symmetry class as H for any A. 

The parametric two-level correlation function is de- 
fined as 

Rp(E+,E-,\) = $>(£- - e-(X))S(E+ - e+ (A))) . 

(2) 

It is mapped onto a dimensionless energy scale, where 
the mean level spacing is rescaled to unity. One has 



A>(r,A) 



lim D 2 R (E+,E-,X) 



(3) 



which solely depends on the difference r = (E + —E )/ D. 
The cross form-factor is obtained as a Fourier transform 



Kp(t, A) = f °° (l - X p (r, A)) e 2 " ltr dr , 



t > 



(4) 
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Time t is measured in units of Heisenberg time tn — 
D^ 1 . Fidelity decay is expressed via the echo operator 

a 

M(t, A) = exp(i2ntH- (\) / D) cxp(-i2ntH + (\)/D) . 

(5) 

Its expectation value with a given state is the fidelity 
amplitude and its average 



(6) 



is a measure for the difference in the two time evolutions 
as a function of A. 

The functions in Eq. ([3]) were calculated exactly with 
the supersymmetry method Q for (3 = 1,2,4. The 
Fourier transforms are (see also [l3j ) 



^i(i,A) 



, " 2t 2 {t-u)(\-t + u) 
du I av— 



max(0,t-l) Jo { v t 2 ) 2 

exp (-2tt 2 A 2 [2u< + t - t 2 + v 2 ] ) 

X y /(u 2 -v 2 ){u 2 + 2u+l-~v s ) ' 

exp (-2^A 2 i 1 + 9 < t - 1 )) / , , , 
' sinh ( 27r 2 A 2 i 2 



2, 2 XH 

Ar 4 (t,A) = i 2 / d w / + ; 



0(t-M ; 



(t 2 - v 2 ) 2 
i - 1 + i) exp (-7r 2 A 2 [f 2 - v 2 + 2tu\) 

y/[(u - l) 2 - V 2 ] [(u+l) 2= ^] 



, (7) 



with Heaviside's ^-function. For A — the cross form- 
factors reduce to the standard form factors Kp(t) 
i.e. Kp(t, 0) = Kp{t). In Fig. [1] we show Kp(t,X) versus 
time t for two values of A. For A = 0.1, the correlations 
vanish as t — > oo. A second peak develops in the GSE 
case for A = 1 at t = 2. The singularity at t — 1 persists. 
For the GOE and for the GUE cases finite peaks appear 
at t = 1 but not at multiples thereof. For all ensembles 
another peak appears for small times t <C 1. Its location 
scales asymptotically with A~ 2 . 

The peak appearing at t = 1 for large A clearly indi- 
cates that here the correlations decay more slowly as a 
function of A than at all other times t. We study this in 
more detail by an asymptotic analysis in A of the exact in- 
tegral expressions ([7]). We calculate the weight Wp{t, A) 

— A 2 ft\-2 Kp(t',A)dt' of the peaks at t = 1 and, for 
the GSE, also at t = 2. In contrast to the peak height the 
weight is well defined for all times for all three ensembles. 
We find 



(8) 



and W4:(2, A) oc A -6 . The weight of the first peak t w 
scales as A~ 2 independently of the ensemble. These 
decays are governed by power laws in A while they are 
exponential for all other times. We shall see below, that 




FIG. 1: (Color online) Cross form-factor (black) and fidelity 
(red) versus time for two different values A. The results for 
the three Gaussian ensembles /3 = 1,2,4 are given as thick 
solid, thin solid and dashed lines, respectively. 



the behavior of the cross form-factor at t = 1 is directly 
related to fidelity revivals, which for the GSE also occur 
at t = 2. 

For the classical ensembles Kp(t), is non-analytic at 
t = 1 |4]. The degree T>(g,x) of non-analyticity of a 
function g(x) at x is defined as the smallest integer T> 
for which the V-th derivative 

i s 

discontinuous at x. 
For the form factor we find V(K i} 1) = 0, V(K 2 , 1) = 1 
and V(Ki, 1) = V(K4, 2) = 3. For typical times we find 
T>(Kp, t) = co, because Kp(t) is analytic. We thus arrive 
at a relation between the asymptotic behavior of Wp(t, A) 
for large perturbation to the degree of non-analyticity of 
Kp(t) which reads 



Wp(t,X) oc \- 2 v(K„t ) + 



t > 



(9) 



We conjecture that this relation also holds for arbitrary 
iM 1,2,4. 

We use the multiply kicked Ising (MKI) spin chain pro- 
posed in 0, [13] to illustrate the revival in the cross form- 
factor. The MKI spin chain is a periodic 1-d array of L 
spins 1/2 with anti-ferromagnetic nearest-neighbor Ising 
interaction of unit strength and periodic boundary condi- 
tions. Each spin receives periodically two different kicks 
of instantaneous magnetic field pulses. The time-reversal 
breaking Floquet operator of the system is Uuki = 
UiU^UjUk , where U\ is the time evolution operator of 

the unkicked spin chain and — exp(— ' &j)> 
(ri = 1, 2) describes each magnetic pulse with a dimen- 
sionless magnetic field b^ n '. <jj are the Pauli operators 
for particle j. The translational symmetry [5j — > 
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FIG. 2: (Color online) The cross correlation function for 
the MKI model with L — 18. Filled/empty triangles cor- 
respond to the real/imaginary part of the cross correlation. 
The statistical error (measured by the imaginary part) is small 
enough to observe clearly the peak. The theoretical expecta- 
tion Eq. (O is plotted as a thick curve. 



which foliates the space in L different symmetry sectors. 
For the choice = (0, 1, 1) and = (1.4,0, 1.4) the 
spectral statistics in most symmetry sectors display ex- 
cellent agreement with the GUE. We introduce an ad- 
ditional magnetic pulse of strength S in z direction as a 

perturbation. We define Us — E/mki ■ ex P ^ J^^a a j ) 
and calculate the cross form-factor of E/mki and Us using 
direct diagonalization, omitting the problematic sectors. 
The perturbation strength A can be calculated from S 
using the correlation functions of the perturbing opera- 
tor Details are given elsewhere. 

In Fig. [5] we compare results of this model with RMT 
results of Eq. ([7]) . We see good agreement with the theo- 
retical result, up to statistical fluctuations, measured by 
the imaginary part. In particular the peak at t = 1 is 
observed. 

We now derive the announced differential relations 
connecting the cross form-factor with fidelity, deferring 
the reader to a follow-up paper for more details. Con- 
sider a general 4-point parametric correlation function 
F a p- ll s(zi,Z2) = {T a (3- n s{zi, z 2 ] H~ , H + )) defined as 



zi - H~ 



a/3 



H+ 



7<5 / 



(io) 

The definition of the angular brackets is now expanded 
to denote either the average over an arbitrary matrix 
ensemble with measure dv(H), or the energy averaging 
over a spectral window of an individual quantum chaotic 
system. We do not require dv{H) to be Gaussian or 
even rotationally invariant. The distribution of V, on the 
other hand, is required to be Gaussian in order to ensure 
the existence of the announced differential relations at 



finite order (see Eq. (fT6|) below). The Fourier transform 
of fidelity amplitude corresponds to F a p.^ a and para- 
metric spectral correlator corresponds to F aa .pp, with 
a summation over double indices. Introducing Hi 2 v ^ a 
fi(Hi : 2 — H ± (X)), Fourier transforming the matrix delta- 
functions and integrating over V, the averages are rewrit- 
ten as 



dA 1 dA 2 dv(H)dH 1 dH 2 



where the symmetry class of the matrices Ai, A2 corre- 
sponds to the symmetry class of H , and multiple factors 
of 2ir are absorbed into the definition of dAi^ 2 . 

The invariance of the flat integration measures dH\ 2 
with respect to independent shifts in H\ and H 2 implies 



tr 



( d 



i) 



\dHi dH 2 



(12) 

The full measure in Eq. (|11[) is also approximately in- 
variant under a simultaneous shift of H\ and H 2 . The 
violation of this symmetry stems from the non-invariance 
of dv(H) under the shifts of H. However, universality im- 
plies that the correlation functions depend on such shifts 
only through the average density of states and level veloc- 
ity variance Q . This dependence is thus manifested only 
on time scales much shorter than tn , which is of interest 
here. In invariant unitary RMT ensembles universality 
under shifts was shown in [l5j . Although not yet proved 
in general, no violations of this universality are known. 
In particular universality follows automatically in models 
which allow for field theoretical representations of corre- 
lation functions [16| . With these caveats we can set 



tr 



d 



8 



dHi dH 2 



f 

J aa; 



00 







(13) 



and combine with Eq. (TT2[) to 

d d \ 1 

4 ( tr 'Qjf'Qjf^ aa '00) = \ tr ( Al - -kz) 2 Faa;Pf3 



(14) 



Using 



tT dH x dH/ aa ^ 



d 2 



dzidz 2 



F a /3;0a, (15) 



4/3 d 



(tr(Ai - A 2 ) 2 T aa -pi3) = -j^Q^Factfp (16) 

and Fourier transforming F, we finally show that 

d ~ Air 2 t 2 
—K {t,\) = — />(£, A), (17) 

under very general assumptions. The averaged fidelity 
amplitude has been calculated in Ref. Q for the GOE 
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and the GUE. For the GSE, 

f +1 Z -1 -!"! (u + t) 2 -l 
U(t,X)= du dv { J ] \v\0(u-l+t) 

J -l Jo \ l ~ v ) 

(i 2 - v 2 + 2tu) exp (-tt 2 A 2 [t 2 - v 2 + 2tu] ) 

X yS[(u-l) 2 -V 2 } [(U+1) 2 ^J • 

A direct comparison of the exact expressions for Kp(t, A) 
obtained in the present contribution and of the ones for 
fp(t, A) in Ref. [3] and in Eq. (| 18[) confirms the validity of 
Eq. (|17|) in the universal RMT regime (although we stress 
that it is valid for any disordered/chaotic model which 
exhibits a separation of scales between the oscillatory lo- 
cal and smooth global behavior of spectral statistics). 

Relation (| 1 T[) allows us to view fidelity revival at 
Heisenberg time tn as being rooted in the algebraic de- 
cay of the cross form-factor. Furthermore, due to the 
established relations, power law decay as a function of 
A must also hold for fidelity at tn and, for the GSE, at 
2tn- This could also have been derived directly from the 
exact equations. 

In Fig.[H we show the fidelity amplitude. Similar to the 
behavior of the cross form- factor, a peak at t — 1 appears 
for all three ensembles Q and for increasing (3 the peaks 
become more and more pronounced. In the GSE case, a 
second peak emerges at t = 2. This peak was not seen in 
the numerics of [3j as it was beyond numerical accuracy. 

Relation (fl~7|) is, essentially, a Ward identity associ- 
ated with the action (fTTj) . It immediately allows to es- 
tablish a connection between fidelity amplitude and the 
Fourier transform of the level velocity correlator C(t, A), 
which is related to if by a Ward identity — 47r 2 i 2 C(i, A) = 
(d 2 /d\ 2 )K(t, A) (see, e.g., 0). As seen from Eq. (TTJ), 
if is a function of A 2 , it follows from (flT]) after a short 
calculation that 

pC{t,\)= (^2 + 4A 2 ^/(i,A) . (19) 

To summarize, we established relations between cross 
form-factor and level velocities on the one hand, and fi- 
delity decay on the other hand. They hold in any sys- 
tem displaying universality of spectral correlations. The 
present formalism can be used to construct a whole fam- 
ily of Ward identities relating apparently unconnected 
correlation functions. One instance are generalizations 
of the 'optical theorem' found in 

[HI 

which relates fi- 
delity amplitude for small perturbations to the spectral 
form factor Kp(t). Further, the results presented here 



do not apply to crossover regimes, where V changes the 
symmetry of H. One such relation was obtained using su- 
persymmetry methods in [n| • A broader set of differen- 
tial relations, generalizing those of (loj . can be obtained 
by utilizing different transformation properties of the ac- 
tion Eq. (|11[) under symmetry-preserving and symmetry- 
violating shifts. Details of these and other hierarchies of 
relations will be presented elsewhere. 

Our findings make it possible to explain features of one 
quantity via the other, i.e. the characteristics of fidelity 
decay in terms of the cross form-factor or vice versa. In 
particular, the revivals of both quantities are linked in 
this way. We studied in detail the decay laws of the cor- 
responding peaks. Further peaks are not possible. The 
very occurrence of the peaks in the cross form-factors is 
neither trivial nor intuitive and will be discussed in else- 
where. 
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